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THROUGH SPREADING 


: 
Over a period of the Angeles County Flood Control District in 
California, ‘under the writer’s : direction, has been engaged in water conservation 


ee the diversion of clear stream flow to off-channel areas, causing the water to a 4 


percolate into t the soil from shallow basins and thereby adding i it to the ground- 


water ‘supply. y. This operation is is known as spreading, the areas upon 


which it is performed are called : spreading grounds? 


In the course of speeding « a s mound will be formed above the initial ground- 


water table. The mound will be two dimensional three dimensional, 
depending on the shape of the ground. ‘The “recharge of 
= through wells is a special, three-dimensional cane. The control of sea 

water intrusion into ground- water basins | is of. particular ‘significance. “The 


configuration of | the mound, its limiting 1 ‘magnitude, its effective 


e, temporary 
Ww water retention both a as to volume and time, its flow net for r potential flow, : and a 
‘the drainage ¢ of the water stored in are analyzed i in this | paper. 


Finally, the theoretical results are compared with results derived from model 
tests on two-dimensional flow. 


The 1 mound is two. denen, if spreading is s done from strip areas and — 
three dimensional if the area approximates : a circle. The height H to which 


- a mound may rise during a spreading operation depends on: (1) Rate (qe) of 


_ spreading—that - is, the rate of infiltration; (2) duration (T) of spreading ; ; 
(8) permeability (K) of the soil; (4) porosity. (pu) of the soil ; (5) vertical dis- 


tance: (D) between the initial pane 4 water surface and the ground surface; 


_ Norte. —W ritten ‘comments are invited for publication; the last discussion should be by 
Asst. Chf. Engr., Angeles Flood Counted Dist., Calif. 
- 2“Correlatin sé Flood Control and Water Supply, Los Angeles Coastal Plain,” by Finley B. Laverty, 2 
Transactions, 4 ASCE, E, Vol. 111, 1946, p. 1127. 


| 
4 
| 
| 
— 
| 
= 
| 
4a. 
= 
| 
= 
= 
+? 


> 


GROUND WATER SPRE EADING 

(6) depth (a,) between initial panera ater surface and impervious stratum; 
(7) slope (¢) of the impervious stratum ; and (8) horizontal distance e (La) to a 
control through which the downstream boundary “conditions are uniquely 

The mound phenomenon must be divided into four distinet p phases. 

- first three occur during its s growth. and the fourth phase duri ing its ;daappearante 


after discontinuation of spreading. 
‘The First Phase.—The first phase consists of a mound wave created w hen 
the inflowing water is being put into storage, necessarily above the initial 
ground-w water surface. ‘For constant inflow this wave must grow 
The growth of 
wave continues until a ‘boundary i is “reached. The in regard to x 
7 is called the “lateral control,” as represented by a a river, a a lake, or a trench. 
The boundary i in regard to y is sealed the ‘ “petentie control. ue It is the surface | 
The Second Phase— —The second phase is transitory. begins when: 
_ (a) The wave front reaches t the lateral control, or (b) the wave crest reaches — 
potential control . Incase 2 (a), the wave length L has reached its maximum 
value. i Discharge now begins and increases until it is equal to the inflow q,. 
In case (6), the height H of the wave } remains constant while the length | LT 
increases and, in the absence of a lateral control, theoretically becomes infinite. 
Hence, all the flow is put into storage at a steadily decreasing rate. - pa 
 ~‘The Third Phase-—In the presence of a lateral control, stationary flow will 
during the third phase. Under these conditions potential flow will 
develop, and the storage in in the mound will reach its maximum for a constant i 
inflow qo It may greatly y exceed available surface storage. 
‘The Fourth Phase.—If the mound has reached pe in the presence 
of a lateral control and if spreading is discontinued suddenly, it will drain 
itself in the course of time. _ If the mound has reached the e potential control 
‘in the absence of a lateral control, and if ‘spreading has subsequently been 
continued until, theoretically, the inflow q Qo has become zero and the wave 
length L has become infinity, then this mound will never drain itself as long 
as there is no extraction from it. This serves to emphasize the fact that the 
‘drainage phenomenon is primarily related to the existence of a lateral control. 


"7 an General Conditions. —The mathematical analysis of the mound phenomenon ~ 


‘Saree 


is based on the Darey 4 of percolation of water through s sand. This law 
is valid as long as the iw is laminar. Laminar flow is governed by aa 
-Reynolds®* ‘number R. M - Muskat? and others have found that, for sand, 


“laminar flow takes place if R 


r _ Furthermore, the idealized soil is assumed to be homogeneous as well as 


q isotropic. Hence, , the p permeability remains constant irrespective of the direc- 


“Les fontaines publiques de la ‘ville de Dijon, ” by Henri Darcy, Dalmont, Paste, 1856. 
4 ‘Recherches physica- -mathematiques sur la théorie du mouvement des eaux courantes,’’ by G. C. F. 

R. de Prony, Mémoires de l'Institute National, Paris, 1804. 

5 “Papers on Mechanical and Physical Subjects,” ai Osborne Reynolds, C J. Clay and Sons, Cam- 


“The Two Manners of Motion of Water,” by Reynolds, Presecdines, Royal Institution « of 

Great Britain, Vol. XI, 1884-1886, p.44. 
Flow of Fluids Through Porous Media,” by M. Muskat, Hill Sek 
Co., Ine. New York, Ast Ed Ed., 1987. 
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conditions 
porosity of the soil and slope of the i imper' vious stratum may be expressed by 


Generally, the storage in the part. of the mound directly beneath the 


spreading ground is small as compared to. ‘the total storage. Therefore, it 


may be neglected except in special cases. 


The slope (2) * the it impervious s stratum is generally small, namely, 


order <0.01. 
PART —THE TWO- —-DIMENSIONAL FLOW OF WATER THROUGH 


TS 
GRANULAR SOM 


1. — ANALYSIS OF Unsteapy Grounp-W ATER 
ai Outstanding among the classical presentations on unsteady ground- wate 


4 flow is the one due to M. J. Boussinesq.' 8 He derived the differential cat 


ay _ Ka 


—w —which expresses this and in which and y are the 


of the free water surface. If the flow becomes steady or stationary, Eq. 1 
. 


rhich i is know n as the Dupuit® equation or theorem. The validity of = 
requires that the vertical of the velocity vectors are small as 


‘on n the Dupuit equation lead to 
Eq. . Lis not integrable by exact calculus methods. Considering the fact 


— 


that, nor mally i in ground- water flow, the slope of the water surface and the — i 


height of the mound relative to the initial ground-water depth | @, are small, - 


in neglecting and in in making y = 


order of magnitude as 


0: 


“Etudes et sur mouvement des eaux Jules Dupuit Carilian- 
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_GROUND- w ATER SPRE ADING | 


negligible 2 with equal justification. Actually terms are of opposite sign 


so long as i i and | SU g are of equal sign. Their difference, then, is of higher order 


smallness. the neglect of both terms is, probably, | a closer approxi-— 
the neglect only. The exception t to this rule and 
nificance will be in its ‘proper place. 1 With this neglect, or for 


= 0, in Eq. 3 there results 


which | is identical i in — the well- known for the 
flow of heat through a nonradiating prismatic bar. Boussinesq ‘referred to 

this analogue three quarters of a century a ago. 
Heat Analogue ‘Solution; Bar of Finite Length” (First, Second, and Fourth 

Phase). —Based on the analogue of the flow of heat through a prismatic, non- 


radiating of length La with a continuous and constant heat source at the 
origin « = 0 (spreading ground) and constant temperature at the opposite | 
- (control), the following solution for the first, second, and fourth 


To typography, let: 
12 
. ‘It may be verified ified through differentiation that | 
> 


GQ et icosBx + Cz 


is is a particular solution of Eq. 4 in which Ci, C3, and are arbitrary con-_ 
stants. It does not satisiy the boundary ‘conditions for the growing 
phenomenon, however. These conditions 


2.y=Oatx=La for all values of ¢ in which 0 < La < ~; and 

Based on the fundamental concept that , y <a, upon which Eq. 4 is | predi- 
cated, t the linear for = “= a 


os 


the equation for the stable at t = ‘The 


is 


 @“Conduetion o of Heat in n Solids,” by and J.C. Yacger, Oxford Unive Press, ‘Oxford, 7 
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“WATER SPREADING» 


‘Hence Eq. | 6 becomes: 


cos 


Eq. 9 still does not t satisfy condition 1 
4 A Fourier series expansion for the equation of the he stable mound lis - 
y - ( i) (La — 2) = ~ 
X cos 


n then becomes 


= 
= 3 7 


The Disappearing Mound Between the Spreading Ground and the 


trol —The boundary conditions now become: — 


all of the conditions. 


\Ka 


(b) : at z = Lg for all values of t; and 


106 indicates thet (a) is satisfied 
8g — K ai) La 


and ‘that (b) (c) are likewise satisfied. ty 
; ‘The General Case of the Unsteady Mound for i ¥ 0- —The e foregoing was 
"predicated on the presence of discontinuity immediately upstream from the 

spreading ground, for example, i in the form of a basin barrier. In the absence ~ 
q of sucha discontinuity the growth of the mound downstream from the spreading - 
7 ground gives 1 rise to a backwater condition ‘upstream from it. H Hence, there 


i 
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— 
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_and the final equati aaa 
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GROUND- WATER SPREADING 


are three distinct zones to be considered as indicated in Fig. 1. ‘The t ae 
configuration of the ‘stable mound will be analyzed in in a subsequent rnd tion. a > 


In the following, the mound surface in zone II is considered parallel to, and 
the the ground- water For a constant 


‘Spreading Ground_ 


Zone III 


Control 


B)= = (He —y) 4% loge 7, 
E> 0: “hal 


7 in which B = width of spreading ground. — 2> 


Ha 


which 0 < 7 < 0.01. 
The stable ‘mound will now be considered approximated by Eqs. 7 and 14, 


mound i in one equation, the origin at oe epreading 


6 
3 
4 
inear equa- : 
1, an approximate, linear eq 
— % 


> 


GROUND-WATER § SPREADID 


which La = L. An equation which 1 satisfies 


For the g1 growing ‘mound two equations are are necessary, for 2 zone III— 


| for I— 

La 


=0. 
16 and 17, with the previous ‘equations defining La, 
stitute a general solution of the shape of the unsteady ‘mound within the limits 
of the initial assumptions and approximations made. _ Fig. 2 shows t the con- 
figurations a at various times ¢ for the growing (dash aren and hers disappearing - 
= 4 ft; qo = = 0.82 
cu ft per hr; K-= = 25 ft 1 per er hr; a = 0.5 ft; and id p= 0. 15 for a | growing: mound 
(dash line in Fig. 2) and 0.06 ler’ a draining ‘meud (solid line in 1 Fig. 2). : These 
values correspond to those i in the comprehensive model tests, except that a #0. 
It is interesting to note a division of flow from the spreading ground between 
zones I and III during part of the growth. . This conforms to the requirement — 
of synchronous growth of the two zones atx =0. The initial ground- water 
is checked during this period. and goes tales storage in gone I. 
This f feature is of particular significance in the control of sea water intrusion. 
It shows that the latter is checked when the infiltrated fresh water reaches 
the saline ground-water table. To maintain control indefinitely the rate of . 
7 infiltration Thust be : such as to cause rise of the m mound to at least mean sea : 


level. The use of recharge wells for ‘this purpose will be ‘in the 


; Eq. 4, which forms the basis for this s solution, holds true strictly f for i= 0 
only, although it should lead to satisfactory for so long as 
: ( uy and 5M are small , of the same order of magnitude and opposite in 


sign. The latter condition is not satisfied during just the 


iq. 4 can be written for the disappearing mound. 
XK 
— 
2e 
-4 
-. 
— 
— 
= 
q 
a 


‘initial per period ol infiltrated flow into zone i. _ However, the influence on the 


‘mound configuration | of this. discrepancy i is o net important in the light of the 

‘mathematical simplification thus achieved. Hence, mathematically, the base 

of the nieaeeeae mound is horizontal. This gives rise to a fictitious outflow 


n Spreading G 


8 

at the upstream e1 end and a discontinuity in slope because of the approximation = 


of the upstream curve by a straight line. To conform to the physical ‘system, — 

the mathematical system ‘must be superposed on the sloping base. “The 
accumulated sum of the ‘mathematical outflows a at the two ends corresponds 

to to o the effective pore \ volume in the mound, drained | ‘ined during time t. t. 


] 
Henee, the out 


he outflow at the control ¢ disappearance i is sa 
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|For I= La Eqs. 16, 17 17, and 1 19 9 express the conditions for the the symmetrical — 


‘True Conricuration OF THE STABLE Mounp FOR 0 


_ After. rotating the z a-axis so as to be parallel to . 


‘The quar iy is neglected i in Eq. 21a. After i introducing Eqs. (21 


Eq. 20 | and 


+ + =2 i| 
Itis approximately true 


when = 0: 
an ye He. 


4 


— (Hat — a%) — iz’ He | 


q 
_ Based on the coo — 
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which is linear, fractional, nonhomogeneous, differential ‘equation. 
stituting andb=),+ N in 28 and solving for M and 
result in the homogeneous equation: __ 
Sy 


Introducing y and db; = 8 wa au + U ds;; integrating and neglecting 


BK 


© 
in whi which (to. simplify 


ry (Aa + 


VBE 
termined by Eq. 12 


12 and shown in Fig.l. 
APPRoxmMATE Equations For Unsteapy Grounp-WaterR 


the Unique Flow FPunction—First Phase (Zone III).—In view of the limi- 


tations imposed upon Eq. 3, it is axiomatic that the respective solutions are :. 
valid only within these limitations, namely, y K a, and — and ¢ is very small. _ 


‘phenomens, "the former quite often is not.  Henee, it to find ¢ 
: unique function g = f(z) for unsteady flow which is in close e agreement with a 
those based on — 3 within the imposed limitations but w which } permits « extrapo- 


g into zone mm; 


must fade away ast —> 


| 
| 
| 
— 
} 7 s 
y= 
lz: 
— The backup curve in zo: . Hence the configura | 
30 
itions: 
lq=q 
q= 0 


— — 


The first three conditions ar are satisfied t by a multitude of | functions, whereas — 
an n exponential function is : required to satisfy condition 4. All four conditions , 


In accordance with Fig. 3: 


Surface 


‘After rotating axes SO SO as as to make the z-axis parallel toi: t: 
aa 
since yi may y be neglected again and 


> 
K Ca + (x’ + L’) ty — vas dx’ = C.. (387) 


on ‘the left side (which & 4 a small fraction of the ar area 


ne 


. — 
q 
— 
4 
- 
Des 
— 
— 
— 
— 

Wety’\, 


=h’? +2 a, + = 


Fora’ =L',y' = 0, and 
| 


in zones II and I. Except for a short initial period, the storage in zone Il 
| be neglected. ¢ follows from the condition of of synchronized rise 
ata’ = 0. . Forti = = 0, and neglecting : zone II, = do. 
Case (a) —At the instant the mound wave _ reaches the lateral control * 


2 = L’, the discharge there is zero, and the configuration of the wave is 
expressed Eq. 40. ™ the other h hand t the configuration of of the: stable mound 


The qu quantity flows into zone III remainder of. Yo storage } 


then” within this period ‘the term eL’ - _ in Eq. 40 


must fade away; h’a must change to Ha (which is accomplished i if (2 ai 
Eq. 41 fades a away); and | q’ must change to Qo. Furthermore, : at any r time t the 
total inflow must equal the total storage plus the total outflow. = =~ 
The unique time function that ‘satisfies these conditions i is: 


and t the: mound equation for the second te is 


in 
in which: 


| 
Case (b) (b).- —It ti -If the mound wave om to the potential control in the absence 
a lateral control, then from Eq. Ai, , it becomes evident that { for the 
_ height h’, to remain constant, the product q’ L' = = Jo -L’ must remain constant. 
Hence, as the wave continues to grow in length, the inflow gq. must diminish 


proportionally. i is strictly. true for t=0 nearly s so for iL’ K2 Qo. 
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Ya _ From a practical point of view this significance of the wave reaching the 
ve ‘surface is is that “waterlogging” of 1 the ground may become critical, not — 


First Phase.—Since the initial ‘ground- water Qo 
changed throughout the wave configuration, it is evident that the water stored © 

in zone III must be caused entirely by the flow q’ and that, in order to satisfy - 
continuity: 


The Approximate Time Requirement.— 


in which: V’; = pore volume; and T’, = time ‘required | to create wave length fo 
-L’ and height h’. The exact integration of Eq. 45 would be a formidable task © 
so long as 7 is a factor in Eq. 40 and q’ is variable. A graphical solution would 
be called for whereby the relation between t, q’, h’, F , and L’ , for. arbitrary slopes i . 
could be established with close approximation. 
Actually the influence ¢ on the final result of 0<i < 0.01 is such that it 
not differ significantly from that due to i= = 0. 
For = 0 Eq. 40 reduces to: 


K 


and Eq. 41 to 


‘The determination of the pore volume of the mound w ave as represented 


th 


by Eq. 46 is facilitated by e expanding it into a binomial —-, converges for 


= 


It that, for a = 0, h? from Eq. a7 i is twice that Eq. 48. 


‘Therefore Eq. ATi is closely by 


tq L ice = 


& 
* — 
— 
| — 
| 
(46) 
> 
4 
— 
an a 
|  Fory = 


- 48 if would to be a function of z for a, = <= 
However, for , Eq. 48 is a approximation if ifé= 1. the . 


Ratio’ 


and upon integration 


¢€ 


A. 
‘The Second Phase, Case (a). —The binomial expansion of Eq. | 43 for 


outflow qatxz= L at any time, ¢, is: 
we 
@= 
and the volume AV, of water discharged during t is: 


ince t is not a fected in n regard to 


from Ba, Eq. 55 


+ 
14 GROUND-WATER SPREADING» 
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_GROUND-WATER SPREADING 
’ ~ which hy» = is constant in order to make integration pene 


Cn i is s determined from Fig. 4 for. The use of this m mean ve value i is seciaaiais 
since the ratio h/ varies between the order of 1/2 and 1/3 for a = 0° and 
ris a very large value | of Qo; and é also applies to the binomial of Eq. 25 for 


=Qandz=L. For ‘normal values of a, the ratio h/H is nearly 
namely 1/3; and ~ 2/3 Kq. 57: 


and the approximate, total accumulated outflow during t the time an 


Vu =H Ay dz.. 
which. Avi is the between y from Eq. 25 and from Eq. 43 fort 
«Using t the binomial. of both equations: 


KG +a) as) / 


tock 


K (h + a) 


~ 


‘The Unique Time Funetion—Fourth F Phase.— 
based on Eq. 4, of the drainage of the mound 1 upon sudden 


‘spreading after stationary condition has been treated in 


to results which are a fair of the drainage phenomenon for 
During the ‘drainage period the flow control is 


=¢ the 
flow through that section will be zero. ~ ‘Hence the at that must b be 


, an — 
.(60) 
2 
— 
— 
— 
— 
i — 
- — 
— 


— 


GROUND-WATER SPREAD 


zero. Making the assumption that the slope varies linearly with 


. 


“which i is a parabola, symmetrical about the y-axis. 


For For one and the same slope (2 ). = --<. -, the mound for stationary 


flow will have a greater volume than will the mound defined by Eq. 66. Booty 


fore, there is a transition period during which the seaaaed remains constant. 7 
From 25 and =0 the volume is: 


4 


whereas, for the ‘mound shape defined Eq. 66, pore volume 


= 
= 


The « change it is 


nd 1 since during this period the outflow i is to be the time is 


- Following ‘this time interval the tate of of outflow varies according to the 


which satisfies the that its tot 0 


, is equal to Ve 


7 
The configuration o of f the ‘stable mound ; for #0 0 was in in Sect. 2. 


Inthe following, po potential flow w will be treated fori = 0. 
Maximum Potential Flow." satisfy a phys ically possible system of 
the flow net between the spreading ground and the lateral control, motion 
“must be irrotational, and there must exist streamlines and a velocity poten- 


tial Furthermore, the condition of must be satisfied, which 
expressed by the Laplace e e quation 


= 

“Hydraulik,” by Ph. Forchheimer, B. G. T ner, Leipsig ‘und Berlin, 1930 Pp. 82-83. 


= ‘Theoretical Soil Mechanics,” by Karl Terzaghi, John Wiley & Sons, Inc., New York, N. Y., 1943, 


and 


— 
| 
a | 
ig 
1 
— 
| 
— 
| 
| 
— = 4 
| 
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ta 
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— 
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— > 


th 
parabola. The ‘equation for for the phreatic line is 


ER 

In order that this sy system shall reprosent f aap under the action of gravity, 

| 


—+y 
—_— 


An analytic comy complex function,” with the added condition that 2% = K 
; along the phreatic line, satisfies the ma mathematical requirements of potential — ; 
flow. _ To satisfy t the theory strictly, the surface of the spreading ground and 
; of the control should be at a constant potential. or condition i is — 
not satisfied, a situation which leads to distortions. 
a. Changing» the constants in the Kozeny™ to pn 
4 
ana 
— w=ot..... 


a potential a stream function; and 


in phreatic line is defined by It is the basic. 


13 “Fluid Mechanics for Hydraulic by Hunter Rouse, Mc Hill Book Co., ‘Ine., 
York, N. Y., 1938, pp. 96-98. 
“Waseerkraft und "by Kozeny, 1931, Vol. 26, 28. 
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0: 


Along ¢ the flor, = =0, ; 80, Eqs. 77 and 73: 4 


"which i is the pre pressure » head 61 on ‘the floor. Generally, y Ke y and 
in which n can be a any number between 1 and Hence: 


Fre. 


Therefore, the velocity for any streamline at -_ point is expressed as: 


— 


tr 


Also, at any point on streamlines 


Knty 


— 
Focus, F for family of conjugate parabolas~ 2k | 
4 


GROUND-WATER SPREA SPREADING 


dy 


which | means 1s that a at the e control, the slope a every streamline is 45°. x. For 


(88) 


‘point. The horizontal transport of the streamline is- 

veosa = —nK 


which i isa maximum when 


Hence, sin a = cos a, which is tow for a =4 ‘80 for the maximum dis- 


charge through the least area the aes must be 45°, which i is the condition | 


This flow ne net represents, in pry — flow far from a — ground > 
and in particular the case of a narrow ‘spreading strip from which the — 7 


merely rains” down upon the mound below. 
‘The case of streamline flow from the spreading ground to the — is 


much more ‘involved. Apparently no rigorous mathematical solution to the he _ 
problem has s so far been found. | ._ Therefore, the problem was | approached as as 
follows: As s indicated i in Fig. 5, _ the complete flow ni net for the basic agra ; 


“alls for the ate through the horizontal seepage e surface between z = 
n 


1. Hence, n reversing and. inverting the parabola, as | shown i 


«Mig. 6, in a double parabola. If a flow g q were infiltrated between F; 
3 and the vertex of the inverted parabola an and dan additional flow q were infiltrated : 


; Hence, at zx = 0; and y = Kn — 
- 
— 
Ss 
tm 


_GROUND-WATER SPREADING DING 


between n the vertex x and of then q = q at the former 
2 2q at the the latter. tter. From 0 and = therefore 
ence, the distance between y=qandy =2¢q¢ along y= = 


+00), 


7-5¢ 


f 


f "Horizontal in 


common to both nets. | Blimi- 


i 


and the double parabola is a streamline, _ Elin 
nation of transforms: this streamline e into the phreatic line. 
a Thus far, each section of the curve has been referred to its own n natural, * -¢ 


coordinate system. _ The equations of the two parabolas | can be referred to 


common system of coordinates i in which the plane of of symmetry will be x = | itd a | : 
the impervious stratum y y = = 0 (Fig. 6). 4 


7 
is; 


= Values of the equipotential ¢ will = at 


7 g 
ty 
11.29 10.97 10.58 |‘10.16*459.73 [49.32 8.51 [%8.13 
ent 
91b) 
* 
4 
) 


y= = 0, the. point of lowest potential the system. So, at y= H= 


: 


a = 


flow net between = and “2 + 
method,’”® that i is, by progressive approximation. — 


The « completed flow net is shown in Fig. 7 using the following oe 


«iL = 46 in. = 3.833 ft; and A = 12 in. = 1.00 ft. Hence from Eq. 91 


Therefore, = ~=>+=0.10 ft =1. 20 


_ 3.833 0.0667 ft. 
‘in.; and ® = K = K 


“Height, in Inches 


Pa 


we 


Ps 


‘system to the final approximation of the relaxation method. - 
‘tial stagnation is” s indicated below a boundary ‘approximated by a 
- ellipse (Fig. 7). ‘The parabolic boundary indicated corresponds to { the one 


_ General Potential Flow. —The conditions for the general case of pc potential 
y VE 1, between q/K and H—are indicated i in He. 6. 7 
‘Three distinet stages of a , must be ¢ considered, namely : Qos 

<a, Qo < D; and (3) D < = H, of which the first case has already 


been “treated. For this case from 91b- and from 


1/2: 


+ q/K,y = 


_, ee umerical Solutions for Tile Drainage of Layered Soils,’’ by James N. Luthin and R. E. Gaskell, 

Transactions, American Geophysical Union, Vol. 31, 1950, p. 595. 
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‘For the second case, , at ;downstreamend: = 


At Ly= therefore: 


‘The for both parabolas is at y= 


yields 


“For the third case, z= L, in Eq. 91a and 


is zero when D = 


‘The flow net for the downstream parabola for = S a, still follows from 
77 and 78, the focus, howe ever, ‘shifting the distance— 


_ beyond the control. At 0 


distance from vertex; hence at x = - 


— 


Condition at the Control.— —Fundamentally the arises because of 
~ change i in pressure along the | seepage surface OA, , Fig. 8, relative to that of the 
basic parabola. _ The latter pressure is shown as a cross-hatched area A. 
‘Under | this pressure condition the flow between the ne surface AB 
If the sand from OA is pressure condition is 
= changed. — If a, = 0, then the pressure at Sect. OA is atmospheric throughout — 7 


and the p pressure gradient between AB and OA is increased. - Based o on hori- 


‘ 
which lies at a 
which hes at a 
— 
— 
= 
: 
| 
4 
4 
® 
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GROUND- -WATER SPREADING 
| = flow paths as a sities yet fair approximation, the | flow g through surface | a 
| OA ~ 1. 50 qo to satisfy the Darcy law. . But since the flow q. is ‘stationary, 
the s seepage surface OA must adjust itself to the new pressure condition and, 7 
therefore, must reduce its height. The critical height 
{ 


b 


A: With the same > deg gree of approximation as before it may be shown that : a 
 tailwater depth = 0.60 hy produces p pressure on Sect. which is similar 
to the pressure, , represented by area in Fig. 8; that generally 


that for this condition OA ~q./K. Finally, it is that for a ground 


* depth a, > a, Sect. OA must be greater than q./K so that a seepage 


= 


surface will always m: maintain itself the raised, line and the 
_ tailwater surface. For this condition the pressure gradient between A A”B 


Values of 2 


in which = = horizontal ength of flow path which ‘not ot differ 


: significantly from that for "area ABO = A; of the basic parabola. . Hence, A 


4 
| — 
— 
7 U9) 
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WATER 


For a, > K the height of then seepage surface rapidly w with 
a,andmaybeexpressed by 


3 —THE THREE- DIMENSIONAL FLOW THROUGH 
_ GRANULAR SOIL 
area by : a “circle, te be treated as three 4 
sional, as shown in ‘Fig. 9. ioe unsteady flow radiates uniformly from ‘the 
‘ “Surface of spreading basin 
II —Using a unique flow. function, similar that used for or the two- 
on leads to | 
+ at, — log. + (x L-R) — R (108) 4 
al — | log. = + (L— +R 


| 


¥ 
i 


g rate per unit area. 


25 
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“Initial Inflow, (%) 


4 100 200 400 500 600 700 


Percentage of Initial Length, 


_ For zone IL and for a, >> H the binomial expansion of Eq. 108 may be used 
2L 


th ‘the ‘volume of wa water stored i in n the 


‘+ +5 ) J. ...(114) 
Recharge Well—First Phase.—If, instead of : a spreadin ng ground, 
ala g 
a well is used in einen of ground water, then R relative 1 to L, and the storage — 


‘int the v well, , may be neglected. .. ‘Then approximatel 


in which Aq is the spreadin Integrating: 

For r=R,s=h and, therefore, from Eqs. 109 and 

: 
| 
as 
which 7”” = total spreading time. Hence: — 


GROUND- WATER SPREADING 


Circular Spreading Ground—Second Phase. —It would not be reasonable to 
consider the analogue of case (a) as previously analyzed, in connection with | 
the three- -dimensional mound because of the improbability of the ‘existence of - 

-acircular control. z However, the analogue of case (b) is all the more important. 

_ When the three- dimensional mound has risen to the potential control, it 
follows from Eq. 109 that for A to remain constant L and gq, must be related - 
at all times: as shown in Fig. 10, curve a. It is based on typical numerical 
values, namely, R = 500 and initial length L= 5,000 ft. 


Recharge’ Well—Second Phase.—As R is as L, the second 
"term under the ‘radical. in 109 reduces to. (loge R/L +0. 72) Ww hich 
se a typical, initial ratio R/L - = 0. 01 leads to the relation shown i in Fig 10 


as curve 


al It i is possible to to gain an insight - into | this complex ease by a assuming’® the 


spreading flow as originating from a source of strength m, discharging vertically 
onto an inclined plane of infinite width and length. If the source is assumed 


at the e origin : and if U = initial uniform ground-water flow w per linear ft, parallel 


to the z-axis, then: en: 


a then the combined flow of the initial uniform m flow, U, and the spreading flow g ~ 


—Uy-n — mtan- 
‘The Sow net based on t these equations may } be drawn by Rankin’ s ‘graphical ; 
‘method! 7 and is shown in Fig. ae boundary of — thound is given b by 


streamline y= 0; so, from Eq. 119 


‘Theoretical Hydrodynamics,” by L. M. Milne-Thomson, ‘Macmillan Co. Ine. New York, N. 
“Elementary Demonstrations of Principles Relating to Stream- lines,” by Ww. Macquorn 
_ Miscellaneous Scientific Papers, Charles Griffin and Co ndon, 1881, p. 522. oe we 
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GROUND-WATER SPREADING» 
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Lo 
Y Substituting for U and m in aad 122 gives the mound width at z = 
‘Thus, the mound the effect of a training: wall and 
to the initial ground-water stream of infinite width, the mound acts like an 


island of finite width and semi-infinite length. Pine 


Streamlines for Kagi 


= 


W= =2ho=. Asymptotic 


Limit of 


for one of the existing the area of which be approx 


mated by. that of. a relatively small circle, namely: R= 
= * = 100 ft per hr; and hr. Hence for i = 


= = 10, 000 ft ft; and fro from m Eq. 123.2, to = - 
In this case, based on the potential funetion the rise of the mound | above the 


normal ground- -water surface is of the order of 35 ft. . This i is ‘its ultimate 


height. 


w 


and that the m S upstream to the stagnation point 
. ong ii 
— 
| 
— 
— 
— 
— 


GROUND-W 


= 0 and L = 5,000 ft, from Eq. 112, 7 = 33.47 ft Bq. 114 
‘Vv’ = 34.23 X 10° ft* = 786 acre ft and tT” = 34.23 hr of spreading. fs 
_ This clearly brings out the magnitude of ‘the mound for i = 0 even during 


‘short spreading operations. often extend over much longer periods. 
PART III. —THE -MODEL TESTS ON TW 0-DIMENSIONAL 


of GROUND-WATER FLOW 
1.—Description oF THE MopEL SeruP 
_ Two diatinet model arrangements were used to test the mechanics of the 
four phases for two-dimensional flow. The The ‘model sand was clean, hard, 
stream-worn material which had passed a n a no. 30 screen and had been r retained 


ona no. . 40 —. grain size, therefore, varied li little from 0. 5 mm. — —_ 


second ‘the outflow hydrograph during the second 
drainage of the mound during the fourth phase. A general view of the model 
setup in connection with this arrangement is shown in Fig. 12, and details of : 


the model proper are » shown i in Fig. 13. - latter consisted of 0: ie ay 


watertight box approximately 12 in. wide, 12 in. and 46 in. long, 


“filled with sand with a grid- lined glass front; 


A ‘metal frame, 4 in. wide at the m end, simulating the spreading 
ground and extending clear across the box; 
Twenty-five piezometers s arranged i in horizontal sets of five, each a at five 
elevations above the bottom and each set extending 2 in., 4 in., 6 in. 8 >. 
and 10 in., respectively, into the sand from the end wall; 
“Eight deep wells, each equipped with threaded brass. rod and steel con- 
tact point, adjustable as to vertical position within 0.005 in., , and connected 
to a direct-current t low voltaye > source; 
- an e. Nine shallow wells behind and parallel to the deep wells, and equipped _ 
with one movable contact point | of the same design as previously described ; 
f. Nine argon bulbs to indicate visually the instant of contact between the a 
"ground -water surface and the point gages; 


@ 


=  * A lateral control comprising a a 40-mesh vertical screen, a tail bucket, — 


‘An automatic comprising x nine active supersensitive relays 

and one spare supersensitive relay equipped with arms and pens for recording | 
_ of the aforementioned well-point contacts in terms of time - latter being 


controlled by: a separate clock, ‘Tecording 1-sec time intervals) 


An automatic stage recorder to ascertain rate af in terms of 


j. A constant- head t. tank with float 1 valve contr ol, fed from a a large ‘primary 


tank and connected to a ‘municipal source of water. supply; 
ke An orifice-equipped manometer to determine the rate of inflow to the 7 


j ‘model. _ The rate of inflow was CO controlled by : a needle valve, and the water 
: 
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GROUND- WATER SPREADING 

was was delivered to the spreading ground by a swing spout. The latter served to 

take volumetric measurements when desired. 


ground were ¢ accomplished by a a trough, c connected = manifold with twenty- 


four small spouts placed at the surface of a }-in. layer of coarse diffuser sand — 


um aeration and uniform distribution of inflow ‘into ‘the spreading 


a] 
- 
= 
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Values of H in Inches 


Wave Height of Mound, in Inches. 


No 
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Edge of Spreading Ground 


= 13.25 Min, Eq. 63 - 
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Surface; 


ymmetry 


~ Water Surface in Weir Box 
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Plane of S 


Values of h in Inches 
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the of the on flow. For this purpose 
deep : and shallow wells were removed, together with the electrical recording 7 


apparatus and the water-stage at the 


— 


TTT TT 
Outflow Hyd h, 


~ Experimental Hydrograph Adjusted 
for Inflow Fluctuations Only 


Outflow, in Cubic Feet per Hour 


5 6 7 8 11 «#12 «13 «214 
Time, in Minutes After Commencement of Spreading 


— 


Outflow, in Cubic Feet 


¢ 


Drainage Hydro h, Eq.71~4 
inag ydrograp Gg. 


‘Time, i in n Minutes After Cessation of Spreading 
ra. 


‘The depth of the sand was increased to 123 in., , the length kept at at 46 in., 
and the s spreading | ground was reduced to 1.20 in. in width, in | line with the | 
in the theoretical study (Fig. 7); eight piezometers were 
twenty-five in n the mas wall, and 


extended in in. into. the sand. 


overlying the test sand. A small U-tube connected to the bottom of the — - : 
spreading-ground frame served to detect air lock, if present. a 
— 
y 
9.901 
0.80 
: — 
0.60 
0.30 
__ Experimental Hydrograph Adjusted for | | 1 {| 
0.10 Weir Forebay Storage and Inflow Fluctuations — a = 
2 17 18 19 20 
— 
oF 
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GROUND- WATER SPREADING 


‘Next, the sand was removed from the surface down to within 0.5 in. of —_ 
theoretical phreatic: in accordance with Fig. 7, near the 


Finally all the sand was parabolic cons 
with the focus of the phreatic line and the vertical end wall as a tan seman 
with: a transition to the horizontal floor near the midpoint. The sand — 


ip 


replaced to 0.5 in. over the theoretical phreatic line to conform to — 


L TEst Runs 
All four phases of the mwend phenomenon were included in the — cnnie. 


_ The rates of flow, actually observed i in the model (which - has the neat width 


3 of 11. 25 in.) are distinguished from the r: rate of flow per linear foot, by Q. 
Rates of flow corrected for a width of 12 in. are shown by g. | The latter are 
based, therefore, » on the same premise as t the theoretical rates of flow. 

For the purpose of recording the first and second phases by the chrono- 
graph, the point gages in | the wells were connected to the positive pole, and the 
pond water was connected to the negative pole of a direct-current circuit. | 

OA comprehensive | test of the first two phases was conducted for a ground- 


‘witha stage a, = 6 i in. ii ‘consisting of ‘many consecutive runs at constant 


Qo =0. 82 cu ft per hr. - During this time all the eight contact points were set 


* “il at progressively higher levels to obtain the configuration of the 1 unstable mound 
_& above the undisturbed ground- water surface in terms s of time. ne. =~ 


Spreading ground box box — 


Height, in inches 


093 


‘Length, in 
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| 


various s initial ground- stages a, Qo W ith the prime 1e objective of 
Go (max) for each stage u under this condition . It was not possible to maintain 


a continuous pressure betw the ground (B/2 = 4. 00 in. 


The recordings \ were by direct, volumetric 
_ After rearranging the model as previously outlined, the test runs on poten- 


With the 12 i in. above the bottom; 7 


(2. With t the | sand surface shaved down to within | 0. 0.5 in, of the theoretical 


— With | the sand surface the same as in condition 2 and in addition | with 


Fluorescein was injected ‘into the spreading basin, and at intervals ot 


: In addition, dye w was inserted adjacent to the , glass is various elevations of he 
length, and streamlines "were drawn from the movement of this dye. 
‘Under this condition | there was no difficulty | in maintaining a continuous 
‘pressure between the ‘spreading ground =1, 20 in.) and the 
contro. 


3.—CoMPARISON Between THEory AND TEST 
Comprehensive Tests.— —Numerical values ‘corresponding to the ‘model tests 
were introduced into the equations. Since the floor « of the model was hori- 


: zontal, only those equations based on i = 0 were considered, “except for the 
- outflow from the mound during its growth and disappearance which was also © 


computed for = 0.01 01 for the sake of comparison. 
~The Reynolds number for the test sand, based on a unit gradient anda 


water temperature of 60° F, was 0.4, Hence, the flow was well 
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GROUND-WATER SPREADING 
= Gre 
|—Water surface=12.8in. 
_ bo» 
Spreading grounds 4 


Top of sand 
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Height, in Inches 
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— Water surface = 12.8 in. 


—s min nin) Travel fringes at time indicated 
ne Pressure (potenial) lines 
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4 indicated by piezometers 
Top of san and by 
° 
A 


Top of sand and atmospheric 


pressure line coincide 
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3 Te Travet tringes at inicated 


= Pressure (potential) lines as indicated by | piezometers 
>> ae 


Height, in Inches 


420 
laminar range. The Permeability coefficient K of the model sand averaged ~ 


25 ft per hr at 60° F . A value of des 0.15 represents average conditions of 


"porosity for the rising stage; and = 0.06 represents average conditions of 
7 porosity for the falling ‘Stage w within the Tange 0 of the mound of the compre- 


Figs. 14 to 17 show the evaluation of Eqs. 7, 10b, 19, 25, 43, 44, 46, 47, 52, 


68 55, 56, 63, 70, and 71 pertaining to all the four phases. — Numerical values: 


used were a, = = 0.5 ft; Go = = 0. 82 cu ft p per hr per unit width; K = 25; sand p 
equals a value between 0. 15 and 0. 06. The rate of propagation of the wave 2 


front during the first phase was gc apenend by a a contact- t-point setting of 0.005 in 


wey 


7 


Coarse sand 


— 
4 — 
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/ | x Pressure (potential) lines as 
10 
58 
=) 

= 

4 || 
— 
ob 2 0 12 
above the initial ‘ground -water surface. The observed average rate over the 
length of the model was approximately 330 ft p per hr. It is closely checked : 
q the theoretical rate (Fig. 14). Based on K= 25 ft per hr and an average 
‘slope of 0.015, he corresponding a average velocity of steady flow would 
0.875 ft per hr or nearly one nine hundredth of the average wave velocity. - 

: _ This explains the : surprisingly short time intervals, often observed in the field, : “1 
between beginning of spreading operations a and the ris rise of the water surface : qv, 
As previously stated, the effect of storage AVm in the mound directly 
beneath the spreading ground was neglected in the derivation of the a : 
for the unsteady mound wave. Therefore, correction was 


accomplished graphically. as shown Fig. 15. For the observed ‘ratio 
0.0328 and = 0.15 this of the theoretical wave height he 
- eliminates the major part of the discrepancy between experimental and theo- — 

: retical values. Since the flow into the mound wave was affected by aati 
in 1 addition to the storage AV, a flat check could not be expected. . Both 


influences” diminish with an increase in h and, therefore with Vie In a strict 


sense, , Fig. 15 is applicable only to the first phase, 


Potential Flow Tests.— —Fig. 18 shows the test result for various” rates of 
bee rved stationary flow Q for condition (1), which was necessarily affected 
— flow. Hence the agg of flow shown at corresponding atmos atmospheric | 
ressure lines included the latter. 
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GROUND- WATER SPRE. ADING- 

-. Fig. 19 shows ‘the reeult for condition 1 for ion 3.52 ft per hr, K = 24 ft 
per or hr, and = 12.80 in. Streamlines, equipotential lines, and fringes 
are indicated. Flow through the “entire prism, except 1 near the control, is 
evident. © Fig. 20 shows the result for condition 2 for g = 2.20 cu ft per hr © 
and K = 22.6 ft | per hr. Hence, q/K = 0.097 ft? as. ‘compared to 0.067 ft?, 7 
the per linear foot being understood i in this discussion of the 
for H = 12 in. on which Fig. 7 is based. | Capillary flow is evident in the small 
wedge of sand next to the spreading g ground which could not be removed for 


Teasons of stability. Partial stagnation in the lower left-hand corner is 


$#$#}B “10 a2 46 


¢ spicuous, and so is concentration of flow immedi: below the phreatic 
= in the upstream half of the model. _ The effective head H averaged | 12.80 
in. and, therefore, was in excess of the theoretical one as shown in ‘Fig. 7. — ‘The _ 
extra head was necessary to assure positive pressure on the floor of the senting > 
ground. The phreatic line maintained itself above the theoretical one as would — 
be expected, ¢ except within the aforementioned wedge. — ° Removal of even the 
smallest amount of sand from the parabolic s surface below the wedge caused 


free water to appear. r. There wa was no difficulty reasonably 
_ Fig. 21 shows the result for condition 3. ~ Replacing t the ‘sand after installa- 
tion of the parabolic floor resulted in a ‘permeability which “was greater. than 
in the preceding test—namely i= 24 ft per hr. - Hence, for the corrected. 
- observed ¢ q = 1.95 ft® per hr q/K was 0.081 ft?. It is interesting to note that 
_ the phreatie line, contrary to. the preceding test, dropped below the surface 
of the sand where it remained parallel ¢ to the surface. | Therefore, the parabolic — 
‘boundary tended to give rise to a similar rate of flow as for condition 2 through 
smaller cross sections. Although tests for an elliptical, lower boundary as 


in ‘Fig. were not it logical to presume that 


In n general, the results the model if 
for * capillarity i is made, fall close to the theoretical ones and, therefore, may be | 
considered a of the theory of the mound phenomenon. 
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‘They writer is indebted to H. E.. Hedger, M. ASCE, chief engineer of the _ 
Angeles County Flood» Control the paper and 


_ offering practical suggestions ; to Maxwell F. Burke, Assoc. M. ASCE, who 


‘supervised the model tests; to Karl F. Smallin and Milan W. Ransom, ‘Assoc. 7 


: M. ASCE, whose mechanical skill and perseverance were instrumental i in the 
' achievement of the comprehensive test results; and to Richard W. Rutowski 
for conducting the tests on potential flow, for checking the mathematical work, 
_ and for his assistance and initiative in the evaluation of equations. _ 
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